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The exponential rate of convergence and the Central Limit Theorem for some Markov 
operators are established. The operators correspond to iterated function systems which, for 
example, may be used to generalize the cell cycle model given by Lasota and Mackey lll2|] . 


I. INTRODUCTION 


We are concerned with Markov operators corresponding to iterated function systems. The 
main goals of the paper are to prove exponential rate of convergence and establish the Central 
Limit Theorem (CLT). It should be indicated that the first result implies the second. The opera¬ 
tors under consideration are more general than those used by Lasota and Mackey in in. The 
authors studied therein some cell cycle model, in which the rate of convergence is already eva¬ 
luated by Wojewodka iliil] . Hille at el. Jl] proposed the generalization of the model and assured 
the existence of a unique invariant distribution in it. We have managed to evaluate the rate of 
convergence, which provides asymptotic stability at once, as well as allows us to show the CLT. 
The results bring some information important from biological point of view. To get more details 
on biological background of the research, see Tyson and Hannsgen 120 ] or Murray and Hunt B- 

In our paper we base on coupling methods introduced by Hairer in |3]. In the same spirit, 
exponential rate of convergence was proven in B for classical iterated function systems (see 
also 10] or B). However, we use coupling methods not only to evaluate the rate of convergence. 
It turns out that properly constructed coupling measure, if combined with the results for statio¬ 
nary ergodic Marokv chains given by Maxwell and Woodroofe llTH . is crucial in the proof of the 
CLT, too. If we have the coupling measure already constructed, the proof of the CLT is brief and 
less technical than typical proofs based on Gordin's martingale approximation. What led us to 
this intriguin g so lution was an unsuccessful attempt to follow the pattern given by Komorowski 
and Walczuk lllll . It is worth mentioning here that an auxiliary model, described by some non- 
homogenous Markov chain, is needed to take adventage of coupling methods. While reading the 
paper, one may see that it is a bright idea to express the Markov operator of interest by means of 
an auxiliary one. 

Similar approach may also help to establish the Law of the Iterated Logarithm (LIT). The proof 
of the LIT is supposed to be provided in a future paper. Some ideas useful for proving it may be 
adapted from Bolt et al. |0]. However, we strongly believe that using an appropriate coupling 
measure will, again, make the proof much easier. 

The organization of the paper goes as follows. Section 2 introduces basic notation and defini¬ 
tions that are needed throughout the paper. Most of them are adapted from H, B,B and B- 
Mathematical derivation of the generalized cell cycle model is provided in Section 3. The main 
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theorems (Theorem [T] and Theorem O are also formulated there. Sections 5-7 are devoted to the 
construction of coupling measure for iterated function systems. Thanks to the results presented 
in Section 8 we are finally able to present the proofs of main theorems. Indeed, the exponential 
rate of convergence is established in Section 9 and the CLT in Section 10. 


II. NOTATION AND BASIC DEFINITIONS 

Let {X, g) be a Polish space. We denote by Bx the family of all Borel subsets of X. Let B{X) 
be the space of all bounded and measurable functions f : X ^ R with the supremum norm and 
write C{X) for its subspace of all bounded and continuous functions with the supremum norm. 

We denote by M{X) the family of all Borel measures on X and by Mgn(27) and Mi{X) its 
subfamilies such that fi{X) < oo and fJ.{X) = 1, respectively. Elements of Mi^{X) which satisfy 
g^{X) < 1 are called sub-probability measures. To simplify notation, we write 

f{x)Kdx) for f : X ^ R, e M{X). 

Jx 

An operator P : Mi^{X) —^ Mfin(A) is called a Markov operator if 

1. P(Ai/ri -|- A 2 /r 2 ) = AiP/ii -|- A 2 T’/i 2 for Ai, A 2 > 0, g,i,g .2 C dTfin(A); 

2. Pfi{X) = pl{X) for // e Mfin(A). 

Markov operator P for which there exists a linear operator U : B{X) ^ B{X) such that 

{Uf, fi) = if, Pg) for / G B{X), /i G Mi^{X) 

is called a regular operator. We say that a regular Markov operator is Feller if U{C{X)) c C{X). 
Every Markov operator P may be extended to the space of signed measures on X denoted by 
Msig{X) = {^1 — /i 2 : /ii, M 2 € For /i G Msig{X), we denote by ||m|| the total variation 

norm of g,, i.e., 

IImII + 

where and come from the Hahn-Jordan decomposition of m (see |@]). In particular, if m is 
non-negative, ||m|| is the total mass of m- For fixed x G X we also consider the space Ml{X) of all 
probability measures with finite first moment, i.e., Mj^(X) = {m G Mi{X) : Q{x,x)fi{dx) < 00}. 

The family is independent of choice of x G X. We call m* C dTfin(X) an invariant measure of P if 
Pm* = M*- For M G ATfin(X), we define the support of m by 

M = {ic G X : g,{B{x, r)) > 0 for all r > 0}, 

where B{x, r) is an open ball in X with center at x G X and radius r > 0. By B{x, r) we denote 
a closed ball with center at x G X and radius r > 0. 

In Msig{X), we introduce the Fortet-Mourier norm 

IlMik = sup|(/,m)|, 

/e£ 

where £ = {/ G C'(X) : |/(x) — /(^ < p (x, y \ /(x)| < 1 for x, y G X}. The space Mi(X) with 
metric ||mi — M2 lie is complete (see la], ll^ or lEh. 

We say that the sequence of Borel measures (Mn)nGV C Mfin{X) converges weakly to the 
measure m G Mfin{X) if lim„^oo(/! Mn) = {f, m) for all / G C'(X). It is known (see Theorem 11.3.3, 
10]) that the following conditions are equivalent 
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• {nn)n£N converges weakly to //, 

• lim„^oo(/, Ain) = (/, lA for all / G £, 

• lim^^oo ll/in - Ailic = 0, 

where {nAn&N C Mi{X) and fi G Mi{X). 


III. MAIN IDEA AND THEOREMS 


Recall that {X, g) is a Polish space and let (0, J", Prob) be a probability space. Fix T < oo. 
We consider a stochastically perturbed dynamical system. The state of Xn, for every n £ N, is 
determined by the formula 


^n+l — SiXfiiin+A' 


We make the following assumptions. 


(I) We consider a sequence {tn)neN of independenf random variables defined on (ff, T, Prob) 
wifh values in [0, T]. Distribution of tn+i conditional on = x is given by 


Prob(An+i < i\^n = x) = j p{x,u)du, 0 < A 

Jo 




( 1 ) 


where p : X x [0, T] [0, oo) is a measurable and non-negative function. In addition, p is 
normalized, i.e., fg p(x, u)du = 1 for x £ X. 

(II) Let S :Xx[0,T] —s-Xbea continuous function which satisfies fhe Lipschifz type inequality 
Q{S{x,t),S{y,t)) < X{x,t)g{x,y) for x,y £ X,t £ [0,T], (2) 

where A : X x [0, T] —[0, oo) is a Borel measurable function such that 


a := sup / X{x,t)p{x,t)dt < 1. 
xgx Jo 


(3) 


(III) supjg[o,T] t),x^ < oo for some x £ X. 

(IV) We assume that p satisfies the Dini condition 

rT 


f 


\p{x,t) -p{y,t)\dt < u}{Q{x,y)) for x,y G X, 


(4) 


where u : R+ —)• i?+, u;(0) = 0, is a non-decreasing and concave function such that 


/ 


0j{t) 


dt < - 1-00 for some a > 0. 


We can easily check that if C < 1, we have 

OO 

ip{t) : = E uj{CAt) < -1-00 for every A > 0 


(5) 


n=l 


and limt_j.o p{t) = 0. 
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(V) Function p is bounded. We set <5 := viiix<EX,t&{o,T]P{x,t), M := sup^.gjjf^^gjo.r] 
require (5 > 0. 

We further assume that, for each A G Bx, 


Prob(xn+i G A) := Hn+i{A) and Pfin = Pn+i, 


where 


PfiiA) = 


lx 


lA{S{x,t))p{x,t)dt fi{dx). 


In JH the proof of asymptotic stability is given for the model, while the exponential rate of co¬ 
nvergence is established thanks to some coupling methods in 12211 ^ . 

Without loss of generality, we may think of (X, g) as a closed subset of some separable Banach 
space P[. Then, trying to describe some intercellular processes more precisely, Hille et al. |@] 
proposed a more general dynamical system 


^n+l — Si^Xnpn+l') T Hn+li 

where {Hn)neN, Hn G Ff, is a family of independent random variables with the same distribution 
given by a measure , which is independent of S{xn, in+i) and its support stays in 5(0, e). 

For this reason, we need an additional assumption 

(VI) Let e* < oo be given. Fix e G [0, e*]. Let be a Borel measure on H such that its support is 
in 5(0, e). For every x G X, we set 


I'xi-) = - x)- ( 6 ) 

We assume that S{xp) + h ^ X for every t G [0, T], x G X and h from the support of 
The Markov chain is given by the transition function II^ : X x Bx —>• [0,1] of the form 

ne(x,X) = ^ p{x,t)v%^^^-^{A)dt. 

Then, we may write the Markov operator P^ : Mi(X) —>• Mi(X) as follows 

Pep{A)= [ Ue{x,A)p.{dx). 

Jx 

The case of deterministic protein production, i.e., when e = 0, fits to the framework presented 
by Lasota and Mackey llill and the results obtained there. 

Hille et al. |@] managed to show the existence of a unique invariant measure in the generalized 
model, described above. However, stability was not proven. We want to focus on evaluating the 
rate of convergence, which additionally provides asymptotic stability in the model and allows us 
to establish the CLT. The proof of the CLT is given in Section 10. 

Theorem 1. Let p G Ml (X). Under assumptions (I)-(VI), there exist C = C{p) > 0 and <7 G [0,1) such 
that 


P^p - p*\\c < Cq^ for n £ N. 


^ In both papers the results are proven for stronger assumptions. 
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Now, assumption (II) is strengthened to the following condition: 

(IF) Let S' : X X [0, r] —X be a continuous function which satisfies 

Q{S{x,t),S{y,t)) < X{x,t)gix,y) for x,y € X,t e [0,T], (7) 

where A : X x [0, T] —[0, oo) is a Borel measurable function such that 

fT 

A := sup / X‘^{x,t)p{x,t)dt < 1. (8) 

xex Jo 

Note that (IF) implies (II), due to the Holder inequality, and we obtain that a < \/A < 1. Assuming 
(IF) instead of (II) allows us to show that p* G Mf{X) := {p G Mi{X) : {x, x)y{dx) < oo}, 

which is essential to establish the CLT in the way presented in this paper. It is proven in Lemma 
[15] that p* is indeed with finite second moment. 

Now, choose an arbitrary function g : X ^ R which is Lipschitz continuous, bounded and 
satisfies {g, yj) = 0. Let {xi)i^N be the Markov chain with transition probability function and 
initial distribution p G (A). For every n G A, put 

u 9{xi) + . . .+g{Xn) 

Vn ■ !— 

y/n 


and let denote its distribution. 

'In 

Theorem 2 . Let g G (A) he with finite second moment and let he the distribution ofrjn, as defined 
above. Assuming that all conditions (I)-(VI) are fulfilled and (II) is additionally strengthened to (IV). Then 
converges weakly to the normal distribution, as n ^ oo. 


IV. AN AUXILIARY MODEL - BASIC ASSUMPTIONS 

Our aim is to prove exponential rate of convergence for the model given in |@], as it is stated 
in Theorem [TJ The idea is to use coupling methods. However, implementing these methods 
directly to the model given above does not give the expected results. Instead, we fix a sequence of 
constants {hn)noN C H, where /i„ G 5 ( 0 , e) for all n e N, and consider a stochastically perturbed 
dynamical system 

Xn+l — {Xnj tn+l) ■ — SiXn., N+l) /in+l; U — 0, 1, 2, . . . 


Note that 


Xn+l 1! ) ^n+1 J ■ 

For abbreviation, we introduce the symbol 

Xn°+i ■= (ThA- ■ ■ Th2 {Th^ (a^o, ti),t2 )...), N+i), 

where the upper index xq indicates the point from which the iteration begins. 
Let us further assume that, for every A G Bx, 


( 9 ) 


Prob(xn+i G A) := Pn+i(A) and Ph — Mn+1/ 
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where, for arbitrary /i G 5(0, e), 

:= 


lx Ljo 


lA{Th{x,t))p{x,t)dt 


p{dx). 


We maintain all previous assumptions (I)-(VI). Now, for every h G 5(0, e), we consider an 
operator 


r/j(x, t) := 5(x, t) + h. 


Note that, as a consequence of assumption (II), is continuous and satisfies the same Lipschitz 
type inequality as operator S satisfies. 

We also set 


c := sup Q 
te[o,r] 



+ £* 


> sup Q 
te[o,r],iGV 


Thi ix,t),x 


Obviously, c is finite, because of assumption (III). 


( 10 ) 


V. MEASURES ON THE PATHSPACE AND COUPLING 


Set X ^ X and {hn)neN C 5, /i„ G 5(0, e) for all n G N. One-dimensional distributions 
n)(i (x, •), n e N, are defined by induction on n 


Il\x,A) = 5M) 


Il\{x,A)= [ lA{Thix,t))p{x,t)dt 
Jo 


( 11 ) 


Ill_,Jx,A)= / Ul{y,A)ni;^ ,^Jx,dy), 

J X 


where A G Bx- We easily obtain two-dimensional and higher-dimensional distributions. If we 

assume that, for x G X, •) is a measure on X", generated by a sequence (Ilj^ 

then 


/ii,... 5/in+l 


ix,AxB)= j Ul^^^{zn,B)Up-^Jx,dz), 
J A 


( 12 ) 


where 2 : = {zi,...,Zn) and A G Bx^, 5 G Bx, is a measure on Note that 

njj^(x, •))■■■ hS^' ■)' given by (fTTI) . are marginal distributions of ' "’"^^(a:, •), for every 

X G X. Finally, we obtain a family (x, •) : x G X} of sub-probability measures on X°°. 

This construction is motivated by 1^. The existence of measures 11)^ (x, •) is established by 

the Kolmogorov theorem. More precisely, for any x G X, there exists some probability space on 
which we can define a stochastic process with distribution (f)^x such that 


ct>^x{B) = Prob(r G B) := 5) for 5 G 5x-. 

Therefore, 11^ (x, •) is the distribution of the non-homogeneous Markov chain on X°° with 

sequence of transition probability functions (njj^,)jg 7 v and (j)^^ = 6x, for x G X. If an initial distri¬ 
bution is given by any p G Mfin(X), not necessarily by 6x, we define 

= [ n^,,^,...(x,5)p(dx) for5G5x^. 

Jx 
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Definitions. Let a family of probability measures ({n]j,(a:, •) : x G X})jg 7 v be given. For every i G N, 
we can set another family of probability measures ((x, y),-) : x, y G X} on such that 

• forAG Bx, 

. Cl{{x,y),XxB) = I[l{y,B) forBGBx, 

where x,y G X. For every i G N, {Clf{x, y), ■) : x,y G X} is called coupling. 


VI. ITERATED FUNCTION SYSTEMS 

We consider a continuous function 5 ; X x [0, T] —)• X and a sequence of continuous mappings 
given by {Th.)i^N with sequence of constants {hi)i^x established. We assume thatp : X x [0, T] —)■ 
[0, oo) is a non-negative and normalized function. For each A G Bx, we build a sequence of 
transition operators, as we did in (fTTI) . 

Let n G N. Note that, for arbitrary A G Bx, FI^^ f^^f,A) : X —)• i? is measurable by definition. 
Furthermore, ') ■ R is a probability measure, for every x G X. Hence, 

is a transition probability function on the n-th marginal. Thanks to these properties (see Section 
1.1, for every n G N and a sequence of constants {hfi^zx fixed, there exists a unique regular 
Markov operator for which 11^^ is a transition probability function, and it is given 

by the formula 

{Ph„...,h„h){A) = [ Ul_,,Jx,A)p{dx), 

Jx 

where A G Bx, p G Mi{X). Moreover, a dual operator : B{X) —^ B{X) to y^^ is 

defined as follows 

m,...,hj){x) = [ f{y)nl_y,^{x,dy). 

Jx 

Remark 4. According to assumptions (II) and (IV), one may check, although through some tedious com¬ 
putations, that, for every n G N and a sequence of constants (hi) i(zx fixed, 

l|nfei,...,h„(x,-) -FIl^,...,hfiyr)\\c < a^g{x,y) + ip(^Q(^x,y)), 

where (p is given by ©. This indicates weak continuity of the map X 3 x ^ 

Noiv, this property, together with the fact that Pf^ is a regular Markov operator, impies that Pf^ 

is even Feller (see Chapter 6, 

However, these estimates do not give us any proper result about stability of the model. That is why we 
still need to use some coupling methods. 

Repeating the construction from the previous section, we obtain Pff y^^ p for p G MfiX). 
Obviously, for every n G N, Pf^ y^^p is the n-th marginal of P^ y^^ p. 

Fix X G X for which assumption (III) holds. We define V : X ^ [0, oo) to be 

V (x) = g{x, x). 

Lemma 5. For every n G N and a sequence of constants {hi )i^N fixed, if p G Xll{X), then ^ 

Ml{X). Moreover, 

{V,Pf^_y^^p)<a^{V,p) + r^c, 


where c does not depend on the sequence {hfii^x- 
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Proof. Recall that a < 1 and c are given by @ i (ITOl) . respectively. The state is of the form 
Following ([2]), we obtain 

= / ••• / Q{Xn,x)p{Xn_i,tn)p{Xn-2,tn-l) ■ ■ ■Pix,ti)dtn. . .dti n{dx) 

Jx ^Jo Jo 

■ QiXn, xf) + Qix^,x) f„)p(x^_2, fn-l) • • • p{x, ti)dtn ■ ■ ■ dtip{dx) 

Jo 


< 


< 


< 



X Jo 

^ I - - 

... X{Xn-l,tn)KXn- 2 ,tn-l) ■ ■ ■ A(x, ti)Q{x, x) + g(Xn, X^-l) + • • • + g(xf,x) 
X Jo Jo 

p(x^_^,tn)p(x^_ 2 , tn-l) ■ ■ ■ P(x, ti)dtn • • • dtip{dx) 



/ a^Q{x, x) + c(a"' + ... + l)p{dx) 

Jx 


< a^{V, p) + 


1 — a’ 

which completes the proof. 


□ 


Fix probability measures p,i^ ^ Ml{X) and Borel sets A,Be Bx- We consider b £ Mfin(TL^) 
such that 

b{A xX)= p{A), b{X xB) = u{B) 
and 6^^ £ Mfin(X^) such that, for every n € N, 

bl_,JA xX) = {Pl_j,^p){A), bl^_^SX X B) = {Pl_,,^v){B). 
Furthermore, we define V : X'^ —)• [0, oo) 

V{x,y) = V{x)+ V{y) forx,y£X. 

Note that, for every n € N, 


{VX,,...,kJ < a{V,blff^,^J+2c < a^{V,b) + r^c 


(13) 


For measures b £ Ml^{X‘^) finite on X'^ and with finite first moment, we define the linear functio¬ 
nal 

fib) = / gix,y)bidx X dy). 

Jx^ 

Following the above definitions, we easily obtain 

m < {v,b). ( 14 ) 


VII. COUPLING FOR ITERATED FUNCTION SYSTEMS 

On we define the transition sub-probability functions such that, for A, B £ Bx, 

fT 

Qhiiix,y),AxB)= mm{p{x,t),piy,t)}6(^rh.(x,t),Th.(y,t))iJ^x B)dt, i £ N, (15) 
Jo 
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and 

Qhu...,hr.ii^^y)^Ax B) = [ Qi^{{u,v),Ax B)Ql-\,^^_^{ix,y),dux dv), hGN. (16) 

Measures generated by the transition functions defined above are, by convention, denoted with 
the same letter. Every time, the context should indicate what we mean. It is easy to check that, for 
every i G N, 

Qhiiix,y),A X X) < [ p{x,t)ST,,.{x,t)iA)dt =[ lAiThi{x,t))pix,t)dt = Ill.{x,A) 

Jo * Jo 

and analogously Q\.{{x, y),X x B) < njj,(?/, B). Similarly, for n G N, 

Ql^_f,Jix,y),Ax X) < Ul^_,Jx,A), Ql^_f^Jix,y), X x B) < Ul_^Jy,B). 

For b G Mfin(X^), let denote the measure 

xB)= [ Ql^^,„^hA{x,y),A X B)b{dx x dy) for A,B G Bx,nG N. (17) 
Jx‘2 

Note that, for every A,Bg Bx and n G N, we obtain 

'>)('4 X B) = x B)b{dx x dy) 

= / B)Q)^ f^J{x,y),du X dv)b{dx X dy) 

Jx^ Jx^ 

QL^Ai^,y),A X B){Ql_f^J){du X dv) = {Qi^^,{Ql,_hx.b)){A X B). 

(18) 


/X2 


Again, following (fTT]l and (fT^ . we are able to construct measures on products and, as a conse¬ 
quence, a measure b on X°°, for every b G Mfin(X^). Now, we check that, for n G N and 

b G Ml^{X% 




(19) 


Let us observe that 


4’{Qhu-,hJ)= / Q{u,y)Qhu...,h„ii^,y)^du X dv)b{dx X dy) 

Jx^ Jx^ 

/ / g{u,v)mm{p{u,t),p{v,t)}d(TH^{u,t),ThXBJ))(^^^ 

Jo J 


ix^ Jx^ . 


dv)b{dx X dy) 


lx'2 7X2 


/ / Q(Th„iu,t),Th„{v,t))p{u,t)dtQl^^ {{x,y),duxdv)bidxxdy) 

7x2 Jo 

_ _ _ _ 

/ Q{u,v)X{u,t)p{u,t)dtQ^ ^ h - {{x , y), du X dv)b{dx X dy) 

Jo 

<a / Q{u,v)Qh~^ _ i{x,y),du x dv)b{dx x dy) < ... < aP'(l){b). 

For every i G N, we can find a measure R\.{{x, y), •) such that the sum of Ql.((x, y), ■) and 
R\^{{x, y), •) gives a new coupling measure Cl^{{x, y), ■). 
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Lemma 6. Fix i G N. There exists the family {Rf{{x, y),-) : x,y G X} of measures on such that we 
can define 

Clf{x,y),-) = Qlfix,y),-) + Rlf{x,y),-) for x,y £ X 

and, moreover, 

(i) the mapping {x, y) i-;- R\f{x, y),A x B) is measurable for every A,B e Bx) 

(ii) measures R\f{x, y), •) are non-negative for x,y e X; 

(in) measures Clf{x,y),-) are probabilistic for every x,y & X and so ((x,y), •) : x,y G X} is 
a transition probability function on X^; 

(iv) for every A,B G Bx and x,y G X, we get Cl,{{x,y),A x X) = 11^. (x, A) and C\.({x,y),X x 
B) = Fiify,B). 

Proof Fix A, B G Bx- Let 
Rhii{x,y),A X B) 

= (1 - Qifix,y),X^))-\ulix,A) - Ql((x,y),Ax X))(nl(y,B) - Ql((x,y),X x B)) 

if Qlf(x,y},X^) < 1 and R\f(x,y),A x B) = 0 if Q\f{x,y),X‘^) = 1. Obviously, the formula 
may be extended to the measure. The mapping has all desirable properties (i) — (iv). □ 

Lemma[6]shows that, for every i G N, we may construct the coupling {C^. ({x, y),-) : x,y G X} 
for {nj^,(x, •) : X G X} such that Q^.((x, y), •) < C^,((x, y), •), whereas measures y), •) 

are non-negative. Following the rules given in ifTT]) . (fT^ . as well as the whole construction from 
Section O we easily obtain the family of probability measures ((x, y), •) : x, y G X} on 

(X 2 )oo marginals 11^ (x, •) and 11“ (y, •). This construction appears in |@]. Note that, 

for every n G N and x,y G X, ^^((x,y),-), constructed as in (ITT]) , is the n-th marginal 

of h 2 Additionally, {(7^^ ^^((x,y),-) : x,y G X} fulfills the role of coupling for 

/i„(X) ■)■ X G X}. Indeed, for A G Bx, 


Cl,...,H„ii^,y),A X Cl((u,v),A X ^^((x,y),du x dv) 

= xdv) = ... =U^^ j^Jx,A) 

and, similarly, C'^^^ ,^^((x,y),X x B) = ^^Jy,B). 

Fix (xo,yo) G X^ and (hn)n<sN C [0,e). The sequence of transition probability func¬ 
tions i^^((x,y),-) : x,y G X}^ defines the non-homogenous Markov chain 'F 

on X^ with starting point (xo,yo), while the sequence of transition probability functions 
(7^1 ^^((x, y, 0), •) : x,y G X, 0 G {0,1}}^ defines the Markov chain ^ on the augmen¬ 


ted space X^ X {0,1} with initial distribution C^{(xo, yo), •) = '5(a;o,yo,i)(')- If V: where 

x,y G X,i G {0,1}, then 

Prob(^n+i G A X 5 X {1} I T'n = (x,y,i),i G {0,1}) = y), A x B), 


Prob('l'„+i G A X B X {0} | d'n 


(x,y,i),i G {0,1}) = R])^^ f,^({x,y),A x B), 
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where A, B ^ Bx- Once again, we refer to (fTTI) . (fT^ and the Kolmogorov theorem to obtain the 
measure ((xq, yo), •) on x {0,1})°° which is associated with the Markov chain 

From now on, we assume that processes 'F and ^ taking values in X2 and X2 X { 0 , 1 }, 
respectively, are defined on (n,F, P). The expected value of measures ((a^o, 2/o), •), 

^S,h 2 ,...((^o,yo), •) is denoted by 


VIII. AUXILIARY THEOREMS 


Recall that a is given by Q. Fix x G (0,1 — o)- Set 

-Rx = {{x,y) G : V{x,y) < x“^2c}, 

where c is given by (ITOl) . Let d : (X^)°° — )• N denote the time of the first visit in K^, i.e. 

d{{xn,yn)neN) = infjn G N : {xn,yn) G K^}. 

As a convention, we put d{{xn,yn)n&N) = oo, if there is no n G X such that {xn,yn) € Xx- 
Theorem 7. For every C G (0,1) there exist positive constants Ci,C 2 such that 

-£'xo,yo[(® + ^) ^ CiV{xQ,yQ)-\-C2- 

Proof. Fix (xo, yo) G X^. Let T' = (x„, yn)n&N be the Markov chain with starting point (xq, yo) and 
sequence of transition probability functions ( {Cl. ((x, y), •) : x, y G X} ) . Let Fn C F, n £ N, 

be the natural filtration in Fl associated with T'. We define 

An = {u: £FI\ T'j = (xj(w), yj(w)) ^ Xx for z = 1,..., nj, n G X. 

Obviously, A„+i c An and An G Fn, for n G X. In consequence of (fT3]) . as well as the definitions 
of An and Xx, the following inequalities are P-a.s. satisfied in fl: 

'^An^XQ^yflY {,Xn+liyn+l)\Fn\ ^ l,4n (^{^(^' 11 , yn) T 2c) < 1,4„ (o + >r)X(x„, y„). 
Accordingly, we obtain 


/ X(Xn,y„)dP< / V{Xn,yn)dP= X[l/(Xn, yn)|-Fn-l]dP 

f A.n JJ A.n—1 

< / [oX(xn-i,yn-i) + 2c]dP < (a + x) / X(xn_i,y„_i)dP. 

J A.n—\ J A.ji — i 


On applying these estimates finitely many times, we obtain 


/ X(xn,y„)dP < (a + x)"- / l/(xi, yi)dP < (a + x)” ^[aX(xo, yo) + 2c]. 

'An ^A\ 


Note that 


P(,4„) < f x( 2 c) ^V{xn,yn)dP < >ir[ 2 c(a + x)] ^(a + >ir)"-[oX(xo, yo) + 2 c]. 

J An 

Set c := x[2c(a + >r)]“^[aX(xo, yo) + 2c]. Then, P(A„) < (o + x)"'c. Fix C G (0,1). Since d takes 
natural values n £ N, we obtain 


OO OO OO 

+ ^r‘^"PiAn) < + x)-f”(a + >r)”c = ^(a + x)(i-f)"c, 

n=l n=l n=l 

which implies convergence of the series. The proof is complete by the definition of c and with 
properly choosen Ci, 6 * 2 . □ 
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For every positive r > 0, we define the set 

Cr = {{x,y) S : Q{x,y) < r}. 

Lemma 8. Fix a G (a, 1). Let Cr he the set defined above and suppose that b G is such that 

supp b C Cr- There exists 7 > 0 such that 

for 6 and M defined in assumption (V) (see Section\TTT^. 

Proof. Recall that is given by (H)). Directly from (fT7|l . (fTHl) and dl^ we obfain 

= / / 2 / X dv) b{dx X dy) 

J o J 0 

= / / lCa-ari^n-:XV}^'^^{p{Xn-l^'tn),P{Xn-l^tn)} ■■■xum{p{x,ti),p{y,ti)}dtn...dti h{dxxdy). 

Jx'2 '-J{0,T)^ 

Note fhaf Icgn, {x^, xV) = 1 if and only if (ti,..., G Tn, where 

Tn := G (0,T)’" : g{x^,xl) < a”r}. 

Set Tn ■= (0, T)^\Tn. Note that, according to assumption (II), we have 

/ Q{Xn,xl)p{Xn_i,tn) ■ ■■p{x,ti)dtn ... dti < a^g{x,y) < a^r 

'^n 

for (x, y) G Cr- Comparing this with the definition of 7^', we obfain 

cTr / p{xn_.i^,tn)...pix,ti)dtn...dti <a^r, 

Jt 

which implies 


/ P{Xn-l-,tn) ■ ■ ■ P{x, ti)dtn . . .dti < — < 1. 

It « 


We then obtain that the integral over Tn is not less than 1 — j > (1 — §)” =: for sufficiently 

big n £ N, which provides, using assumption (V), that \Tn\ > ■ Finally, 


{Ql,...,hMCTr) > / 5^\Tn\b{dx,dy)>d^ 


7 


J \ M 

If we sef 7 := 6M~^'y, the proof is complefe. 

Theorem 9. For every x G (0,1 — a), there exists no £ N such that 


□ 


ll< 3 ^,h 2 ,...((a^,y)G)ll > for{x,y) £ 


where 7 > 0 is given in Lemma^ 
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Proof. Note that, for every real numbers tt, u G i?, there is a general rule: min{u, + 7 ;|—u > 0. 

Hence, for every (x, y) G and i £ N, we obtain 


mm{p{x,t),p{y,t)} + |p(x,f) - p{y,t)\ -p{x,t) 


(if > 0 


and therefore, due to (ITSll . 


fT 

IIQL((®’y)’')ll + \p{x,t) - p{y,t)\dt 


> 1 . 


For every b G Mfin(X^), due to the Dini condition (see assumption (IV)) and the Jensen inequ¬ 
ality, we get 

\\QhM= [ Qhi{{x,y),X^)b{dxxdy)= j \\Ql{{x,y),-)\\b{dxxdy) 

JX'2 7X2 

> ||5|| - f Lo{Q{x,y))b{dx x dy) > ||6|| -uj{4>{b)). 

7x2 


Then, by (fTHI) . 

'J -A 

>\\b\\-u{f{Ql^b))-...-u{f{Ql^_j,h)). 

Following (frUjl and recalling that oj is non-decreasing, we obtain 

n 

> \\b\\ - 

i=l 


See dSI) to recall the definition of p. Thanks to assumption (IV), we know that limi_j.o ip{t) = 0. 
Hence, we may choose r > 0 such that if g{x,y) < r and therefore a~^(p{b) < ra“^||6||, then 

Er=i^K“V(&)) < p{a~^f{b)) < ^ll^ll. 

If supp b C Cr, then we obtain 

a PO) 

Fix X G (0,1 — «)• It is clear that C C^-i 2 c- H we define uq := min{n G N : 
a’^{x)-^2c < r}, then C„no^-i 2 c C Cr- Remembering that y), •) = 

{Qhr.+i,...,hr,+mQhi,...,hJ{ix, y), •) and using the Markov property, we obtain 


Then, according to (f20]l and Lemma [H we obtain 

iir, 

■' h\ ,... j/lnQ 




WQZ . h.Ai^^y)r)\c. 


2 ^)’ ^ Qhi,...,hnfj(A^y)^Cano:K-i2c) ^ 

2 “ 2 “ ~ 


for {x, y) G K^. This finishes the proof. 


□ 
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Definition 10. Coupling time r : (X^ x {0,1})°° —)• N is defined as follozvs 

r{{xn,yn, On)n€N) = infjn G N : 9k = l for k>n}. 

As a convention, we put T{{xn,yn, &n)n&N) = OO/ if there is non G N such that 0^ = 1 for every k > n. 
Theorem 11. There exist q G (0,1) and Ca > 0 such that 

< C^il + V{x,y)) for {x,y) G 

Proof. Fix X G (0,1 — a) and (x, y) G To simplify notation, we write a = {a + x )“2 . Let d be 
the random moment of the first visit in K^. Suppose that 

d\ — d, dn+i — dn d o 

where n G N and are shift operators on (X^ x {0,1})°°, i.e. Tn{{xk,yk:dk)keN) = 
{xk+n,yk+n,0k+n)k€N- Theorem [ 7 ] implies that every d„ is (7^ ^^ 2 ,...((®> d)> finished. The 
strong Markov property shows that 


Ex,y[a'^ o TrfjXrfJ = for uGN, 


where E^,^ denotes the n-algebra on (X^ x { 0 , 1 }) generated by dn and 'F = {xn,yn)neN is the 
non-homogenous Markov chain with sequence of transition probability functions ({C^ ((x, y), •) : 
x,y G X})jg 7 v By Theorem[7]and the definition of K^, we obtain 



( 21 ) 


We define T{{xn,yn,On)neN) = inf{n G N : {xn,yn) G Ok = I ior k > n} and a = inf{n G 
N : T = dn}. By Theorem [Oj there is no G X such that 



( 22 ) 


Letp > 1 . By the Holder inequality, (ITT!) and (I22|) . we obtain 



k=l 


k=l 


< [C'iH(x,y) + C 2 ]pr]~p 


k=l 



k=l 2 


1 


For p sufficiently large and q = a p, we get 


Ex,y[q '^] = Xa,,j^[ap] < (1 + l/(x,y))C3 


for some C 3 . Since r < f, we finish the proof. 


□ 
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Lemma 12. Let f £ C. Then, there exist q G (0,1) and Cs > 0 such that 

[ \fiu) - {(x,y),-)){du X dv) < q^C^{l + V{x,y)) for every x,y £ X,n£N, 

Jx^ 

where n* : x { 0 , 1 })°° — x { 0 , 1 } are the projections on the n-th component and : X^ x 

{ 0 , 1 } —)• X^ is the projection on X^. 

Proof. For n £ N we define sets 

A. ={t£{X^x{0,l}r: r(t)<^}, 

B. ={t£{X^x{0,l}r-- r{t)>^}. 

Note that n = 0 and An u Bn = (X^ x {0,1})°°, so, for n £ N, we have 


Hence, 


/OOO 

^hiM 


{{x,y), 


/Soo 

^hiM 


{{x, y), OU™ + C^i,h2,...iix, y), Olfin . 


[ |/(u) x dv) 

Jx^ ^ 

Jx^ ^ 

< [ Q{u,v)iU*^2KCr,,h2,...iix,y)r)\An){duxdv) + 2C^^^^^^^^^^^^ 
J 

Note that, by iterative application of (fT9l) . we obtain 



TT* /OOO 

,h2 


,...{{x,y)r)\An){du,dv) 


= mhKicjfj,h2,...ii^,y)r)\An)) 


Then, it follows from (fT^ and (fTT|l that 


n+1 I - 


4>{U*^2lll^^ {C^M.-.i^x, y), OU. )) < aL VJ y(x, y) + 


2c 

1 — a 


We obtain 


[ 

JX 2 


\f{u) - f{v)\{u*j,2Kc^,,h2...ii^^y)^ •))(d« X dv) 


< 0^2 


LfJ 


aL 2 iY(^x,y) + 


2 c 1 


1 — 0 - 


+ y), Bn). 


It follows from Theorem [TT] and the Chebyshev inequality that 


n. 


= C^,,h2,...ii^^y)dr > tt}) = <^“,h2,...((a;,y),{g ^ >q 2}) 


< 


^x,y \q 


< q2C^{l + V{x,y)), 


for some q £ (0,1) and C 3 > 0. Finally, 

/ l/(^^) - /(^^)l(n^ 2 n;C'^;j 2 ((a:,?/),-))(do x dv) < oLfJ 04(1 + H(a:,?/)) + 2 ^ 1 ( 73(1 + H(x,?/)), 

Jx^ 

where = max{a 2 ,(l — o)“^2c}. Setting q := max{a 2 ,g 2 } and Cs := (74 + 2 C 3 , gives our 
claim. □ 
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Remark 13. If g : X ^ R is an arbitrary bounded and Lipschitz function with constant Lg, then, there 
are q £ (0,1) and C 5 > 0, exactly the same as in Lemma^^for which we obtain 

[ \giu) - x dv) < Gq'^C^il + V{x, y)) for every x, y e X,nGN, 

Jx^ 

where G := max{Lg, sup^-gj^ |5((x)|}. 

Theorem 14. There exist q E (0,1) and Cs > 0 such that 


f&c 

= sup 
/6£ 


\\^h^,...,hjx,-)-Ul^ j^^{y,-)\\c < q'^G 5 {l + V{x,y)) for x,y & X and neN. 

Proof The theorem is a consequence of Lemma ll^ It is enough to observe that 

\K^,...,hAx, •) - Oik = sup [ fiz)(Ul_f,^ix, •) - Ku...T^{y, ■)){dz) 

fee Jx 

[ (/ki)-/k2))(n3,2n;(7;^,;,^,,„((x,y),-))(d^i xdz2) 

Jx^ 

Hence, using the argument of Lemma IT^ we obtain 

\\Ul_^,Jx,-)-Ul_^,Jy,-)\\c<q^G,il + V{x,^^^ 

which finishes the proof. □ 


IX. EXPONENTIAL RATE OE CONVERGENCE - PROOE OP THEOREM 1 
Note that we may write 

PeK-)= [ [ ■■■[ ^hi,...,hjx,-)u‘^{dhi)...u‘^{dhn)p{dx). 

Jx JB{0,e) JB{0,e) 

Comparing this approach with Remark 1 and Lemma 1, we see that is Feller and, for every 
n e N,it satisfies the following property 

{V,Pfy)<a^{V,p) + -^. (23) 

1 — a 

We present the proof of Theorem 1 below. 

Proof of Theorem 1. Let pi, /i2 G ^i^st want to evaluate \\Pfpi — Pfp 2 \\c- Let f G C. We 

obtain 


\{f,Pfpi-Pfy2)\ 

■ [ [ fiz)'^hi,...,hj^^dz)u‘^{dhi)...u%dhn)pi{dx) 

3,e) JB(0,e) Jx 


IX JB{0,i 

j 

'B{0,e) 


IJ 

Jx Jb\ 


lB{0,e) Jx 


fiz)'^hi,...,hr,iy^dz)u‘^{dhi).. .v^{dhn)p2{dy) 


lx 


L 



X JB(0,e) JB{0,e) JX 


f{z)lll^^ , hjx,dz)u‘^{dhi) . ..u‘^{dhn)pi{dx) P2{dy) 



X JB{0,e) JB{0,e) JX 


fi^)^hi,...,hjy^ dz)v^{dhi)... U^dhn)p 2 {dy) pi{dx] 
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Now, following the result from Theorem [141 

•) - K,,...,hSy .Oik < + V{x,y)) 


where q and C 5 are independent of choice of /ii, /12 ,..we obtain the inequality 

|(/,P,>1-P,>2)| 



{dhi)... {dhn) yLi{dx)yL 2 {dy) 



B(0,£) 


0 - ^hi,...,hSy^ ■)\\cv"{dhi)... v^{dhn) Hi{dx)fi 2 {dy) 


<q^cJ [ [[ 

Jx Jx '-Jb 

= / / (1 + 

Jx Jx 


^ JB(0,e) JB{0,e) 

(1 + i/(x,y))/ii((ix)/i2(d2/) 


(1 + V{x,y))i^^{dhi)... v^{dhn) Hi{dx)y. 2 {dy) 


where measures /xi, /i2 G 

Now, set := ^ e Ml{X) and fi 2 ■= G Ml{X), for arbitrary m £ N. Note that it follows 
form Lemma 1 that P^y. is with finite first moment if /i G We obtain 



for some constant Cq. Hence, {P^y)neN is a Cauchy sequence in (Mi (X), || Ok)- It is then proven, 
because of completeness of the space, that {P^y)n£N converges in (Mi(X), || • H^). Put y^.{y) = 
lim„_^oo Ped-- As mentioned before, we know that P^ is a Feller operator and this impies that the 
measure y*{y) is invariant. Then, for yi,y 2 G Ml{X) and every e > 0, we have 


\\y*{yi) - y*{y2)\\c < \\y*{yi) - -P^Filk + ll-ffFi - -PrF2lk + lk*(F2) - tT^2lk < e 


for sufficiently large n £ N. Hence, we have the invariant measure /i* := y*{y) which is unique 
in MI{X). We should make it clear that p* G Ml{X). Note that we can take a non-decreasing 
sequence {Vk)k£N such that Vk{y) = m.\n.{k,V{y)}, for every k £ N and y £ X. Fix x £ X. From 
the first part of the proof, we know that (/, P^ 6 x) converges to (/, y^) for every f £ C, which 
means, by the Aleksandrov theorem (see Theorem 11.3.3 in |@]), that P^Sx converges weakly to 
/i*, since both measures are probabilistic. Hence, for all k £ N,Vk £ C{X) and we obtain 



Note that, according to (l23l) . for every n £ N, 

{Vk, P^5x) = a^{Vk, 5x) + (1 - a)-k < a”W(x) + (1 - a)-k 


{Vk,y,) = lim {Vk,P^6x) < (1 - a)-^a 


and, additionally. 
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so the sequence {{Vk, n*))k^N is bounded. Because {Vk)k&N is non-negative and non-decreasing, 
we may use the Monotone Convergence Theorem to obtain 

/ y{y)lJ‘*{dy) = liTo. / Vk{y)y*{dy). 

Jx Jx 

Then, V is integrable with respect to p*, so p* is with finite first moment. 

Keeping in mind that V(x, y) = V(x) + V(y), the exponential rate of convergence to the unique 
invariant measure p* G Ml(X) derives from the following estimates 


ipru- 




l-c ^ 



^”(75(1 + V{x,y))y^{dy)ix{dx) < 


X JX 


where C := 6*5(1 -|- V{x,y))y^,{dy)ij{dx) < 00 for p G Ml{X). Finally, since C is dependant 

only on p, the proof is complete. 

□ 


X. CENTRAL LIMIT THEOREM - PROOF OF THEOREM 2 


Let us first make the follownig observation. 

Lemma 15. if p G Mf(X) is with finite second moment, then {V'^^Pfiy) < 00 and therefore p* ;= 
/i*(p) = lim^^oo Peh has finite second moment. 

Proof. Let p G Mf (A). Fix x ^ X, n > 1. Recall that A < 1/2 and c are given by ((H]) and dTOl) , 
respectively. Moreover, is of the form (0. Reasoning as in Lemma 1, we obtain 




/ / ... g{x^,x)p{xi^_i,tn) ■ ■ ■p{x,ti)dtn .. .dti pi{dx) 

Jx ^Jo Jo 

/ f-f 
if-f 


2g^{x^, xfif) P 2g^{x^, x) i„). ..p{x, ti)dtn ■ ■ ■ dtip{dx) 


Xfixn, in) • • • A^(x, x) + 2gfix, x^) + 4:g‘^{xf, xf) , x^) 


~x\ 


p{Xn-l,tn) ■ ■ ■ Pix, ti)dtn ■ ■ ■ dtip{dx) 

< p) -h 22 c^( 1 + 2A ... -h 2”A*") 


< 2A^{Vfifi)+4: 


1 - 2A’ 


Estimates are independent of choice of sequence {h^ )nGN and therefore Pfp) < 2A"(K^, p) + 
4c^(l — 2A)“^. We take a non-decreasing sequence {Vf)k^N such that Vf{y) = min{A;, V‘^{y)}, for 
every k £ N and y £ X. We know that Pfp converges weakly to p*. Hence, for all k £ N, 
Vf £ C{X) and 

lim {Vf,Pfp) = {Vlpf) < 4c2(1 - 2A)-\ 


so the sequence ((V/^, p*))fcev is bounded. Because {Vf)k<^N is non-negative and non-decreasing, 
we may use the Monotone Convergence Theorem to obtain 

= lim 

fc—fCXD 


so, indeed, p* is with finite second moment. 


□ 
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Let rjn and ^rjn be as in Section I///I In particular, 77 * and 77 * are defined for the Markov chains 
with the same transition probability function IIj and initial distributions p* and 5x, respectively. 
Further, let : X —)• i? be a bounded and Lipschitz continuous function, with constant Lg, which 
satisfies ( 5 , p*) = 0 . 

Central Limit Theorems for ergodic stationary Markov chains have already been proven in 
many papers. See, for example. Theorem 1 and the subsequent Corollary 1 in Q by Maxwell 
and Woodroofe. The following lemma implies that assumptions of Theorem 1 and Corollary 1 
(Hi) are satisfied. 

Lemma 16. Let g : X ^ Rbe a bounded and Lipschitz continuous function with constant Lg. Additio¬ 
nally, = 0. Then, 


n=l 



/1 j. 2 

p*{dx) 

k=0 


1/2 


< 00 


Proof. Note that, by Lemma IT^ and Remark IT^ 


(24) 


n—1 


n—1 


J2{9,Pedx) = [i9,Ps"dx) - {g,p.) 
k=0 k=0 

n—1 .. „ 

X] / / 9{z){Til{x,-)-lfl{y,-)){dz) p^{dy) 


k=0 ' 
n—1 

E 

k=0 ■ 
n—1 


X 


iaizi) - g{z 2 ))iU*^ 2 niC^^ ,^^ {{x,y),-)){dzi x dz 2 ) p*{dy) 


<Y,Gq^C5 {1 + V{x,y))p,{dy). 

k=0 


Then, for every x £ X, n £ N, 


n 1 '\ 'P' r 

V(5,pX)<G'C 5-^ / (l + C(x,y))/i,(dy)<C'9(l + C(x)), 

where Cg := GC^{1 — g)“^(l + fx P{y)T*{^y))- Keeping in mind that p* is with finite second 
moment, we obtain that (pTl) is not bigger than 


^n-3/2[C|(l + 2C + C^7z*)]i/2 <00 

n=l 


and the proof is complete. □ 

Hence, by applying Corollary 1, we obtain that <^ 77 * converges to the normal distribution in 
Levy metri^ as n —)• 00 , which equivalently means that the distributions converge weakly to each 
other (see (j] for proofs). 

Now, the idea of the proof is based on the following remark. 

Remark 17. Note that, to complete the proof of Theorem^ it is enough to establish that ^rjn converges 
weakly to <^ 77 *, as tt, —)• 00 . Equivalently, it is enough to show that lim„_j.oo W^Vn — 4 > 77 *||£ = 0, since 
weak convergence is metrised by the Fourtet-Mourier norm. 











20 


Proof of Theorem 2. Set x,y £ X and choose arbitrary f £ C. Suppose that we know that the 
following convergence is satisfied, as n —)■ oo, 


IR 


f{u)^Vnidu) - / f{v)^hnidv) 


'R 


(25) 


(26) 


Then, by the Dominated Convergence Theorem, we obtain 

[ fiu)<^r]nidu) - [ fiv)<l>7]*{dv) 

Jr Jr 

< [ [ [ fi'<J‘)‘^Vnidu)- [ f{v)^r]y{dv) n{dx)n^{dy)0, 

Jx Jx Jr Jr 

as n —^ oo. Note that, by Theorem 11.3.3 in |@], (f26]l implies that ^rjn converges weakly to ^rj^f 
as n — )• oo, which, according to Remark[T71 completes the proof of the CLT in the model. Now, it 
remains to show (1^. Note that 


IR 


f{u)^hn(.du) - / f{v)^r]y{dv) 

Jr 


f 




g{ui) + . . • + 


^ Ttg’'''’”(x, dui X ... X dUn) 

f + .»(,,*. x...xd„„) 

Jx" ^ y/n j 

(27) 


where Iif"''^{x, •) = e) • • • /s(o e) ^hi f^^idhi)... v^{dhn) is a measure on X"-. We may 

write 


/ / 

Jx^ Jx^ 


f 


g{ui) + ... + g{un)\ _ j-Cgivi) + ... + g{vn) 


n 


n 


(^> dui X ... X dun)Ti]^f;\^ (y, dvix ...X dvn) 


h\ ,... ,/ln 


< 


/ k(- 


g{ui) + ... + g{Un)\ _ J.fgivi) + ... + givn) 


n 


n 


(28) 


((x,y), ■)]{dui X ... X dun x dvi x ... x dvn), 


where n| ^ : (X^ x {0,l})-^(X2x{0,l})- are the projections on the first n components and 
n^ 2 n : (X^ X {0, 1})" —X^” is the projection on X^”. Since / is Lipschitz with constant L/, we 
may further estimate (f28]l from above 


J- 


L 


\g{ui) - g{vi)\ + ... + \g{un) - givn)\\ {'n*x2r.'n.i,...,nC^i,h2-((^^ y)^ ■)j ^^dui X (iui)r=i) 
X] / 2 l9(«i) -5'(^^i)l(n;^2n*C^;,^ ((x,y),-))(d^ii X dvi). 




Now, for every 1 < f < ra, we refer to Lemma [12] and Remark IT^ to obseve that ([28]) is not bigger 
than 


LfG 


2=1 


^ fC^{l + V{x, y)) = n-" 2 LfGC^q^-^{l + V{x, y)). 


Note that the expression above is independent of choice of sequence {hn)n&N and, thanks to this, is 
also the upper bound of (f27|l . We go with n to infinity and obtain ([2^ . The proof is complete. □ 
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